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Abstract 

We study the observability for a heavy Majorana neutrino N along with a new charged gauge 
boson W at the LHC. We emphasize the complementarity of these two particles in their production 
and decay to unambiguously determine their properties. We show that the Majorana nature of N 
can be verified by the lepton-number violating like-sign dilepton process, and by polar and azimuthal 
angular distributions. Lhe chirality of the W coupling to leptons and to quarks can be determined 
by a polar angle distribution in the reconstructed frame and an azimuthal angle distribution. 



I Introduction 



Neutrino experiments, over the past decade, have shown undeniably that neutrinos are massive 
and have large mixing angles [1]. In the Standard Model (SM) of particle physics, neutrino masses 
can be accommodated by a non-renormalizable dimension-5 operator containing left-handed (L.H.) 
neutrinos, vj_, [2]. Such an operator can be generated at low energy by including heavy right- 
handed (R.H.) neutrinos, vr. However, the R.H. neutrinos are gauge singlets and so Majorana 
mass terms should also be present without violating any gauge symmetry. The consequences of 
massive Majorana neutrinos are well-known [3-5], and have been incorporated into many models, 
such as left-right symmetric theories [6]; super symmetric (SUSY) SO (10) grand unified theories 
(GUTs) [7] and other GUTs [8]; R-parity violating SUSY [9]; and extra dimensions [10]. A recent 
review of TeV scale neutrino mass models can be found in Ref. [11]. 

Many of the aforementioned models contain an extended gauge group or Keluza-Klein (KK) 
excitations of SM gauge bosons. We refer to additional vector bosons charged under the U(1)em 
gauge group collectively as "W"' . If the masses of the W and the lightest heavy neutrino mass 
eigenstate, N, are both on the order of a few TeV, then they can be produced in tandem at 
the Large Hadron Collider (LHC). As first observed by Ref. [12], a W' with mass greater than a 
Majorana neutrino's mass allows the possibility of observing the spectacular lepton number (L) 
violating process 



If a W is discovered at the LHC [13], it is obviously imperative to measure its chiral coupling 
to fermions. In a previous work [14], three of the present authors proposed measuring the W' chiral 
couplings to quarks by studying the process 



It was found that the couplings could be establish as being purely left- or purely right-handed by 
analyzing the polar angle of the charged lepton in the top's rest-frame with respect to the top's 
direction of motion in the partonic center of momentum (cm.) frame. 

We now extend this prior analysis into the leptonic sector via the L-violating cascade decay 
of Eq. (1). More specifically, by reconstructing the polar angle of the lepton originating from the 
neutrino decay in the neutrino rest-frame and with respect to the direction of motion of the neutrino 
in the partonic cm. frame, it can be uniquely determined if the W coupling to leptons is purely 
left-handed, purely right-handed, or a mixture of the two. We show that the distribution of the 
angle made between iV's production plane and its sequential decay plane is sensitive to the W 
chiral coupling with the initial-state quarks but independent of the W coupling to leptons. These 
results are demonstrated through a combination of analytical calculations and event simulations, 
assuming nominal LHC parameters. 
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Majorana neutrinos can decay into either leptons or antileptons, and so W' and N may also 
contribute to the L-conserving collider signature 



For completeness, we have analyzed the polar angular distributions of the unlike-sign process and 
comment on the important differences between the L-conserving and L-violating cases. 

This paper is structured as follows: First, in section II, we present our notation for the W 
couplings to SM particles and neutrino mass eigenstates, and list current constraints on both W h s 
and iV's. In section III, we discuss the production and decay of W"s and iV's at the LHC. The 
like-sign lepton signature, pp — > I .d^jj, its reconstruction, and suppressed background are fully 
analyzed in section IV. In V, we propose methods to measure independently the chiral couplings of 
the W to leptons and to the initial-state quarks. Finally, in section VI, we provide a few comments 
on the contribution of W and N to the L-conserving process pp — > W' — > regarding the 

difference between the Majorana and Dirac neutrinos. We conclude and summarize our results in 
section VII. Two appendices are additionally included. The first addresses neutrino mass mixing 
in the context of W 1 couplings, and the second presents a derivation of the matrix element and 
angular distributions for our like- and unlike-sign dilepton signals. 

II Theoretical Framework and Current Constraints 

There are many beyond the SM (BSM) theories containing additional vector bosons that couple to 
SM fermions, for example: left-right symmetric theories [15] with a new SU{2)r symmetry and an 
associated W' R ; Little Higgs models with enlarged gauge symmetries [16]; extra dimensional theories 
with KK excitations [17-19]. Heavy Majorana neutrinos in BSM models are just as common, 
including Refs. [6-10], and models with TeV scale masses [20-23]. 

In this analysis, we assume the existence of a new heavy electrically charged vector boson, 
W !± with mass Myy, and a right-handed neutrino, Nr. We denote the corresponding heavy 
neutrino mass eigenstate as N with mass m^. We stipulate that M\y is of the order of a few TeV 
and My/ 1 > mjy so as the W' —> N£ decay is kinematically accessible by the LHC, but do not 
otherwise tailor to a specific theory. Regarding the parameterization of mixing between neutrino 
mass eigenstates with SM flavor eigenstates, we adopt the notation of Ref. [24], and extend it to 
include coupling to a model-independent W in Appendix A. This parameterization is accomplished 
with a minimum amount of parameters. 
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A W Chiral Coupling to Fermions 

The model-independent Lagrangian that governs the interaction between SM quarks and a new, 
massive, electrically charged vector boson, W', is given by 

£ = -"7f E K + ^ V § KM '^ W* + 9 q L P L ] d* + h.c, (4) 

where iij (dj) denotes the Dirac spinor of an up- (down-) type quark with flavor i (j); y CKM ' 
parameterizes the mixing between flavors i and j for the new charged current interactions; g RL is 
the W"s universal coupling strength to right-(left-)handed quarks; and Pr : l = \ (1 ± 7s) denotes 
the R, L-handed chiral projection operator. 

We parameterize the new boson's coupling to charged leptons with flavor I and neutral leptons 
with mass m m (for the three light states) or win (for the heavy state) in the following way: 

3 



71=1 

3 



,m=l 



>fp L r + h.c. (5) 



9r G?i) 1S ^ e coupling strength to R.H. (L.H.) leptons; Xt m (Ug m ) parameterizes the mixing 
between light neutrino mass eigenstates and R.H. (L.H.) interactions; and Y^v (Yin) parameterizes 
the mixing between the heavy neutrino mass eigenstate and R.H. (L.H.) interactions. Lastly, ip c = 

— T 

Cip denotes the charge conjugate of the field ip, with C being the charge conjugate operator, and the 
chiral states satisfy Pl{iP c ) = (P R ip) c . In Appendix A, our choice of parameterization is discussed in 
detail. Prom a viewpoint of the model construction as discussed in Refs. [1,12,24], one may expect 
that UU\ yyt ~ 0(1) and W\ XX^ ~ 0{m m /m N ). Since we prefer a model-independent 
approach, we will not follow rigorously the above argument and will take the parameters as 

UU\ YY ] ~ 0(1), and VV\ II 1 ~ 0(1O~ 3 ), (6) 

which is guided by the current constraints as presented later in this section. 

In Eq. (5), the W' is allowed to have both independent right-handed (g q R ) and left-handed (g q £ £ ) 
couplings. Subsequently, the pure gauge states W' R and W' L are special cases of W when 

gtf^O and gf = 0, (7) 

and 

g^ = and gf^O, (8) 

respectively. Additionally, the SM W coupling to leptons can be recovered from Eq. (5) by setting 

g R = 0, and g{ = g. (9) 
Here, g is the usual SM SU(2)x coupling constant. 
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B Current Constraints on W 

We list only the most stringent, most relevant constraints to our analysis here and refer the reader 
to Ref. [25, 26] for a more complete review. 

• Bounds from CMS: The CMS Experiment has searched for Wr and heavy N, where 
Mw R > tin, with the i^i^jj collider signature [27], assuming gR = g. With 5.0 fb _1 of 
7 TeV and 3.7 fb" 1 of 8 TeV pp collisions, the present mass bounds for W' R and TV" are 

M Wr > 2.9 TeV (m N « 0.8 TeV) and m N > 1.9 TeV (M Wr » 2.4 TeV.) (10) 

The search for the sequential SM W', W' SSM , decaying into a charged SM lepton plus fa, 
with g' = g, has also been performed. With 3.7 fb _1 of 8 TeV pp collisions [28], the present 
mass bound is 

M WsSM > 2.85 TeV. (11) 

• Bounds from ATLAS: The ATLAS Experiment has also searched for Wr and heavy N, 
under the same stipulations as the CMS Experiment [29]. With 2.1 fb _1 of 7 TeV pp collisions, 
the present mass bounds for W' R and iV are 

M Wr > 2.5 TeV (m N « 0.8 TeV) and m N > 1.6 TeV (M Wr ps 1.8 TeV.) (12) 

• Global Fit Analysis: The effects of a generic Z' boson on EW precision observables place 
bounds [30] of 

Mz'/gz' > 2.7 - 6.7 TeV. (13) 
For Z' and W bosons originating from the same broken symmetry, we expect similar con- 
straints on M-w< Iqw' since 

M w , ~ M z > x O(l). (14) 

• Bounds on Wl — Wr Mixing: Non-leptonic Kaon decays [31] and universality in Weak 
decays [32] constrain Wl — Wr mixing. The present bound for the L-R mixing angle £ [6] is 

|C| < 1 ~ 4 x 10" 3 . (15) 

C Current Constraints on N 

More complete lists of constraints on low and high mass neutrinos, respectively, are available [24,25]. 

• Bounds from Qv(3(5: For uin S> 1 GeV, a lack of evidence for neutrinoless double beta decay 
bounds the mixing between heavy neutrino states and the electron-flavor state at [33] 

Y^\X^L <5xlO- 5 TeV- 1 , (16) 



Ill 



where the sum is over all heavy Majorana neutrinos. 
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• Bounds from EW Precision Data: A TeV scale singlet neutrino mixing with the SM 
flavor states is constrained [34] by 

iVeTvl 2 ,)^! 2 < 0.003 and |Ktv| 2 < 0.006 (17) 

III W and N Production and Decay at the LHC 

For the remainder of this analysis, we consider for our various benchmark calculations only the 
pure gauge states W' R and W' L , respectively given by Eq. (7) and Eq. (8), and with SM coupling 
strength 

9%l = 9- (18) 
More general results can be obtained by simple scaling. Unless explicitly stated otherwise, we take 

M w , = 3 TeV, miv = 500 GeV, \v CKM \ 2 = 1, \V CKM '\ 2 = 1, (19) 

and use the CTEQ6L1 parton distribution functions (pdfs) [35] for all hadronic-level cross section 
calculations. 

Regarding our choice of neutrino mixing parameters, for mixing between L.H. gauge states and 
light mass eigenstates, we use the PMNS parameterization with mixing angles taken from Ref. [25], 
which includes recent measurements of #13, and take 8cp,cti,ct2 = 0. The bounds from Oz^/3/3 
decay are quite severe and discourage collider searches for L— violation in the electronic channel. 
However, neutrino mixing between the mu- or tau-flavor state and lightest heavy mass eigenstate 
can still be considerably larger in L.H. interactions. Therefore, we use 

\V eN \ 2 = 2.5 x 10~ 5 , [V^jvl 2 = 1 x 10" 3 , and |Ktv| 2 = 1 x 10~ 3 . (20) 

These numerical values are in line with Eqs. (16), (17), and (19); and furthermore, mimic the 
observed fi — r symmetry seen in mixing between flavor states and light mass eigenstates. Where 
necessary, for mixing between R.H gauge states and light mass eigenstates, we apply the unitarity 
condition 

3 3 
^|X, m | 2 = l- ^|C/ fm | 2 , for t = e,n,r. (21) 

771=1 771=1 

For mixing between R.H. gauge states and the lightest, heavy mass eigenstate, we apply Eq. (6) 
and take 

\Y tN \ 2 = l, for £ = e,(i,T. (22) 
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Figure 1: (a) The total decay width for W' R (solid) and W' L (dash); (b) the branching ratio of 
W' RL —t N£ + , with subsequent W' R —t Nfi + (dot) and W' L —t N^i + (dash-dot) ratios; and the 
production cross sections at the (c) 8 and (d) 14 TeV LHC of W' R (solid), W' L (dash), W' R -)• Nt + 
(dot), and W' L N£ + (dash-dot). 



A W Production and Decay 

Under our parameterization, the partial widths for W' decaying into a pairs of quarks are 

nw^tb) = ^ KM '\\ 9 i 2 +9 R 2 ) h 



48vr 



1 - xi 



1 + -xi i 
2 " 



(23) 



where Xi = rrii/Mw', and the factors of three represent color multiplicity. Likewise, the partial 
widths of the W decaying to leptons are 



F(W' -> iv m ) 



g R 2 \X em \ 2 + g e L 2 \U £r . 



Y(W'^IN) = (g R 2 \Y m \ 2 + g[ 2 \V m \ 



M w , 
48vr 
\ M w 



48. V 1 "^] ( 1+ 2^ 
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(24) 
(25) 
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Figure 2: As a function of heavy neutrino mass, (a) the total width of TV" superimposed with the 
N — > IW^ partial widths, and (b) the N — > l^W T , N — > I branching ratios for longitudinal 
(A = 0) and transverse (A = T) W polarizations. 



Summing over the partial widths, the full widths are found to be 



M W f 
32vr 

M w , 
32vr 



A + {l-xl)\2 + x 2 t ) + {l-xl)\2 + x 2 N )\j^\Y lN \ 2 + 2 - £ \X lm \ 



4 + (l 



„2\2 



(2 + x 2 t ) + (l-x%)\2 + x%) 1 -Y;\VtN\ 2 +'' 

£=e 



m=l,i=e 
3,r 

3 E 

m=l,£=e 



£rn I 



(26) 



(27) 



As a function of M^/, Fig. 1 shows (a) the total W decay width; (b) the branding ratio (BR) 
of W — > N£, for £ = e, fi, r, defined as the ratio of the partial width to the total W width, T' w : 

T(W IN) 



BR{W' £N) 



w 



(28) 



and the production cross sections for the pure gauge eigenstates W' R L along with pp — > W'~^ L — > 
N£ + in (c) 8 TeV and (d) 14 TeV pp collisions. 

The production cross section of the W' and its subsequent decay to iV is calculated in the 
usual fashion [36]. The treatment of our full 2—7-4 process, on the otherhand, is addressed 
in Appendix B. Since the u-quark is more prevalent in the proton than the (i-quark, and since 
the dominate subprocess of W' + (W'~) production at the LHC is ud — > W' + (du — > W'~), the 
production cross section of W' + is greater than the W'~ cross section. In a similar vein, the mixing 
between L.H. interaction states and heavy neutrino mass eigenstates is suppressed by | V^jv 1 2 ~ 
O(10~ 3 ), whereas the mixing between R.H. interaction states and heavy neutrino mass eigenstates 
is proportional to |Y£tv| 2 ~ Consequently, the W' L — > N£ branching ratio, and hence the 

pp — > W' L — > N£ cross section, is roughly three orders of magnitude smaller than the W' R rates. 
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(a) (b) 

Figure 3: The partonic level process for a heavy W' + production and decay to like sign leptons in 
hadronic collisions. 

B Heavy Neutrino Decay 

A heavy neutrino with mass of a few hundred GeV or more can decay through on-shell SM gauge 
and Higgs bosons. The partial widths of the lightest heavy neutrino are 

2 

r(iV^^W^) =T T = l-\V m \ 2 m N {l-y 2 w f 

T(N^u e Z) =T Z = —^\V m \ 2 m%(l-y 2 ) 2 (l + 2y 2 ) 

T{N^u,H) = T H = -J^\V m \ 2 m%(l-y 2 H ) 2 (29) 

where Wo,t are longitudinally and transversely polarized W's, respectively, and yi = Mj/mjv- The 
decays of the heavy neutrino through a W are not kinematically accessible. The total width is 

T 

T N = Y^ (2(T + r T ) + T z + T H ) (30) 

where the factor of two in front of I^t is from the sum over positively and negatively charged 
leptons. 

Figure 2(a) shows the total decay width (solid) and the partial decay widths to positively 
charged lepton (dashed) normalized to the sum over the mixing matrices. For this plot the mass 
of the SM Higgs boson is set to 125 GeV. The normalized width grows dramatically with mass 
due to decays into longitudinally polarized VF's and Z's and the Higgs boson. Although the width 
appears to be large at high neutrino mass, for mixing angles on the order of a percent or less the 
width is still narrow. 

Also of interest is the branching ratio (BR) of heavy neutrinos into charged leptons: 

BR(N-^£ ± W =F ) = ELe£o+£r) (31) 

TTot 
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Figure 2(b) shows the total BR of the heavy neutrino into positively charged leptons (solid) and 
individually the BR into longitudinally (dashed) and transversely (dotted) polarize W's as a func- 
tion of neutrino mass. The BR's into negatively charged leptons are the same. As the mass of the 
neutrino increases the Z and Higgs decay channels open, hence the branching ratio into charged 
leptons decreases. Since Tq grows more quickly with neutrino mass than IV, for mjy 3> Myy the 
total BR converges to the BR into longitudinally polarized W's. Also, at high neutrino masses 

r « I> « T z (32) 

Hence the total width is 4Fo and, from Eq. (31), the branching ratio into a positively charged 
leptons is approximately 0.25. This is a manifestation of the Goldstone Equivalence Theorem when 
taking win and V^v as independent parameters. 

IV Like-Sign Dilepton Signature 

A distinctive feature of Majorana neutrinos is that they facilitate L-violating processes, and to 
study this behavior at the LHC we consider the L-violating cascade 

<PA) d( PB ) -> W'+ L (q) -> £+( Pl ) N(p N ) -> €+(pi) el(p 2 ) q( P3 ) q'(p 4 ). (33) 

The two diagrams that contribute to this process are shown in Fig. 3. Figure 4 shows the total 
cross section for the like-sign dimuon process as a function of mjy. In it, the solid line denotes the 
pure W' R gauge state while the dashed line represents the pure W' L state. Since the W' R — > Nfi 
branching ratio is larger than W' L — > Nfx ratio, the cross section for W' R is systematically larger 
than for W' L . Additionally, as the neutrino mass approaches the W mass the cross section drops 
precipitously due to phase space suppression. 

In principle, the conjugate process, ud — > W ~, should also possible at the LHC. However, it will 
possess a much smaller production rate because the ud initial-state has a smaller parton luminosity 
than ud. Despite this, all reconstruction methods and observables discussed below are applicable 
to both processes. 

A Event Selection 

Our signal consists strictly of two positively charged leptons and two jets, a fact that allows for 
considerable background suppression. In simulating this like-sign leptons plus dijet signal, to make 
our analysis more realistic, we smear the lepton and jet energies to emulate real detector resolution 
effects. These effects are assumed to be Gaussian and parameterized by 



E JE 



(34) 



10 



Iff 



0.5 



1 1.5 

m N (TeV) 

(a) 



0.5 



1 1.5 
m N (TeV) 

(b) 



Figure 4: Total cross section of pp — > W' — > fi + fi + W~ times W~ — > qq' branching ratio versus 
heavy neutrino mass at (a) 8 and (b) 14 TeV. Solid (dashed) line corresponds to W' R (W' L ) gauge 
state. 
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8 TeV 


14 TeV 


W' L 


w R 


W' L 


W' R 


Reco. without Cuts or Smearing 


4.6 x 10~ 5 


0.046 


9.3 x 10" 4 


0.95 


+ Smearing + Fiducial + Kinematics (Eq. (35)) 


4.0 x 1(T 5 


0.035 


8.2 x 10" 4 


0.71 


+ Isolation (Eq. (37)) 


2.1 x 10" 5 


0.027 


3.2 x 10" 4 


0.50 


+E? + TTljj Requirements (Eq. (43)) 


i.7 x i(r 5 


0.023 


2.6 x 10" 4 


0.42 


+ Mass Req. (Eq. (44)) 


7.2 x 10" 6 


0.012 


2.0 x 10" 4 


0.35 


cr(All Cuts)/cr(Smearing + Fid. + Kin.) 


18% 


35% 


25% 


49% 



Table 1 : Cross section for pp 



L,R 



H + fi + qq' after consecutive cuts for 8 and 14 TeV LHC. 



where a(E)/E is the energy resolution, a is a sampling term, b is a constant term, © represents 
addition in quadrature, and all energies are measured in GeV. For leptons we take a = 5% and 
b = 0.55%, and for jets we take a = 100% and b = 5% [37]. 

After smearing, we define our candidate event as two positively charged leptons and two jets 
passing the following basic kinematic and fiducial cuts on the transverse momentum, p?, and 
pseudorapidity, r\: 

p j T > 30 GeV, p e T > 20 GeV, rjj < 3.0, i] £ < 2.5. (35) 

Table 1, lists the cross sections for Eq. (33) assuming the pure W RL gauge states at the 8 and 14 
TeV LHC without smearing or acceptance cuts (row 1), and with smearing plus acceptance cuts 
from Eq. (35) (row 2). 

The goal of this analysis is to unambiguously determine the properties of W and N . To do 
so, our candidate leptons and jets must be well-defined and well-separated, that latter of which is 
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measured by 

ARij = ^(A^ + iA^y, (36) 

where Afaj and Ai]ij are the difference in the azimuthal angles and rapidities, respectively, of 
particles i and j. Subsequently, we apply isolation cuts on our candidate objects: 

ARfj n > 0.4, ARjj > 0.3 (37) 

for all lepton and jet combinations, where ARfj m is defined as 

Ai?r n = min ARff. (38) 

lJ i=W',N lJ 

In Eq. (38), the subscript i = W',N on denotes the identified parent particle of l{. The effects 
of the isolation cuts applied at both the 8 and 14 TeV LHC are shown in the third row of Table 1. 
To understand the origin of these precise numbers and parent-particle identification, we digress to 
succinctly connect properties of our chiral Lagrangian to the final-state kinematical distributions. 

B Characteristics of Kinematical Distributions 

Our signal suffers from a very evident ambiguity: either lepton can originate from the neutrino 
decay. The origin of each lepton must thus be determined in order to fully reconstruct an event. 
As noted in section III, the width of N is narrow. Consequently, there is a very small probability 
for the phase space of each diagram in Fig. 3 to overlap, meaning that the interference of the 
two diagrams is negligible. In fact, in the W' R case, the interference is exactly zero because the 
charged lepton from the N decay is left-handed while the charged lepton from the W' R is right- 
handed. Furthermore, since the two diagrams add incoherently, it is reasonable to expect that 
only one diagram contributes at a time. Intuitively, this means that only one of the two following 
momentum combinations will closely reconstruct the heavy neutrino mass: 

mjjj = (pi + p 3 + Pa) 2 or mljj = (p 2 + P3 + Pa) 2 , (39) 

where p% and p$ are the momenta of our final-state jets. 

After calculating both permutations of mjv (Fig. 5), the appearance of the N mass peak is 

,Re 
N 

T- t" V~i f^i T" /~\ fir n I /-\ nnUr t» /~\ r\ *^\t t/~\ "i * n /w _ 

N 



stark. Using the central value of the mass peak, m^ eco ', we identify the charged lepton from the N 
decay as the charged lepton from our candidate event that most closely recovers m^ eco -, i.e., 



Am min = min \rriijj - m§ eco -\, (40) 

where niyj for i = 1, 2 is defined by Eq. (39). 

Independent of reconstructing N, the charged lepton associated with the W decay can be 
identified by analyzing the transverse momentum, px, distributions of our final-state objects. In 
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Figure 5: Invariant mass distribution of m^jj and m£ 2 jj, where i% for i = 1,2 could originate from 
either the W or N. The cuts in Eqs. (35) and (37) as well as the energy smearing have been 
applied. 

Fig. 6, the px distributions of the charged leptons (a,b) and jets (c,d) for the W' L (a,c) and W' R 
(b,d) gauge states. As expected, the lepton identified as originating from the W has a Jacobian 
peak around Mw/2 for both the W' L and W' R cases. To understand the other distributions, we 
consider spin correlations. 

Figure 7 shows the spin correlations of the process in Eq. (33) with the single arrowed lines 
representing momentum direction and double arrowed lines spin. The direction z is defined as 
the direction of motion of the neutrino in the W rest-frame. Each column indicates the spin 
and momentum of the particles in their parents rest-frame with the first column in the neutrino 
rest-frame. Note that for the W' R iyV' L ) the heavy neutrino is in a mostly right- (left-) handed 
helicity state. Hence, for the W' R (W L ) the neutrino spin points with (against) the z direction. The 
decays of the neutrino through longitudinal W are shown in Fig. 7(a) and 7(b) for W' L and W' R , 
respectively, and the decays through a transversely polarized W are shown in Fig. 7(c) for W' L and 



As shown in Fig. 2, 500 GeV neutrino preferentially decays into longitudinally polarized Ws. 
We therefore focus on that case. For the W' R , the lepton from the heavy neutrino decay moves 
preferentially along the z direction. Hence, the boost into the partonic cm. frame will be along the 
charged lepton's momentum. In the W' L case, the charged lepton moves in negative z direction and 
the boost into the partonic cm. frame is against the lepton's momentum. Therefore, the lepton 
from the heavy neutrino decay is harder in the W R case than in the W' L case. The contribution from 
decay into transversely polarized Ws is in the opposite direction. However, as noted previously, 
this contribution is smaller than the decays into longitudinally polarized Ws. Similar arguments 
can be made to explain that the two jets are softer in the W R case than in the W' L case. 



Fig. 7(d) for W' E 
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Figure 6: Transverse momentum distributions for (a,b) the lepton identified as originating from 
the W (dashed) and neutrino (solid), and (c,d) the hardest (dashed) and softest (solid) jets in 
pp — > W — > production. The W' L case is represented in (a,c) and the W' R case in (b,d). 

The energy smearing has been applied. 



As previously stated, identifying well-separated objects in our event is paramount to measuring 
our observables. For 14 TeV LHC collisions, Fig. 8 shows (a) the separation between the two jets, 
Ai?jj, and (b) the minimum separation between the leptons identifed as originating from the heavy 
neutrino and W and the two jets defined by 



min AR 

fc=l,2 



(41) 



where i = W' for the lepton coming from the W' and i = N for the lepton coming from the 
neutrino decay. The solid lines are for W' R and the dashed lines for W' L . The ARjj distributions 
peak at low values for both the left- and right-handed cases. This is due to the W from the heavy 
neutrino decay being highly boosted and its decay products therefore collimated. Also, as can be 
seen from Fig. 6, in the W' R case the lepton from the neutrino decay is harder and hence the SM 
W softer than in the W' L case. Since the SM W is less boosted in the right-handed case, the jets 



14 



W'r:N 



< 



A 



A 



(a) 



( 



A 
I \ 



W' L :N / 



4 



< 



(c) 





z 




TV 










=> 



Wn 



A 



A 



(b) 



A 



W' R :N 



IP 



7' 







==> = 





> 



(d) 



Figure 7: Helicity and spin correlations in the chains Nlr — > £ + W~ — > £ + qq' from decay in 
(a), (c); and from W' R decay in (b), (d). Figures (a) and (b) are for longitudinally polarized SM 
W's, and Figs, (c) and (d) are for transversely polarized SM Ws. The decay goes from left to 
right as labeled by the particle names. The momenta (single arrow lines) and spins (double arrow 
lines) are in the parent rest-frame in the direction of the heavy neutrino's motion (z) in the W 
rest-frame. 



are less collimated and the ARjj distribution has a longer tail for W' R than for W' L . Also, since 
the neutrino is highly boosted, its decay products are expected to land opposite in the transverse 
plane from the lepton from W' decay. Hence, Ai?™™^ peaks near tt for both the the left-handed 
and right-handed case. Finally, Ai?™° is peaked near 2m^ /En ~ 0.7 for both the W' L and W' R 
cases. The AR distributions at the 8 TeV LHC are peaked at similar values, but are more narrow 
than the 14 TeV distributions. Based on these arguments, we define the isolation cuts given by 
Eq. (37). 

The isolation cuts more severely affect the W' L cross section since the ARjj distribution is 
strongly peaked at low values for W' L . As the mass of the W increases, the SM W from the heavy 
neutrino decay will become more boosted. Hence, the two jets will become more collimated and 
the effects of the isolation cuts will be even more significant. Since we will only be interested in 
the angular distributions of the lepton, it is possible to relax the ARjj cut and look for one or two 
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(a) (b) 

Figure 8: (a) ARjj distribution and (b) AR 1 ™ distributions for both the lepton identified as 
originating from N and W . The solid lines are for the W' R case and dashed lines W' L . Energy 
smearing has been applied. 

jets with two like sign leptons. Also, the separation between the lepton and jets from the heavy 
neutrino decay depend on the ratio of itin/Mw'- As mjq/My^i increases (decreases) the lepton and 
jets become more (less) well separated. 

C Background Reduction and Statistical Significance 

The SM background for our £ + i + jj signature has been thoroughly studied for the 14 TeV LHC by 
Ref. [24]. The largest background to our process was found to be from ti events with the cascade 
decays, 

t^W + b^£ + u m b, i^W"b^W~cv m £ + , (42) 

and was also found to be greatly suppressed by the lepton isolation cuts in Eq. (37) . The background 
can be further suppressed by noting that leptonic ti events contain a final state light neutrino and 
therefore a considerable amount of missing transverse energy, $t- This is in direct comparison 
with our signal where all the fx is due to detector resolution effects. The for our like-sign 
leptons + dijet events is shown in Fig. 9 for both the right- (solid) and left-handed (dashed) W 
cases. Furthermore, the two jets in our process originate from a SM W whereas the jets in the top 
background do not. Hence Ifix and dijet invariant mass, rrijj, cuts are also applied: 

$ T < 30 GeV, 60 GeV < m# < 100 GeV. (43) 

The effect of these cuts on the signal rate are seen in the fourth line of Table 1. 

Having obtained a measurement of uin from Eq. (39) and M\y> from the W"s Jacobian peak, 
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has been applied. 

if desired, invariant mass cuts on rri£ N jj and s can be imposed to further isolate the signal: 

\ m lNjj ~ m N\ < 0.1 rriN and \s — M\yi\ < 0.1 Myyi. (44) 

The effects of these cuts are shown in the fifth line of Table 1. 

As y/s increases from 8 TeV to 14 TeV, the percentage of events passing the selection cuts also 
increases. See the final line of Table 1. In particular, we note that relatively fewer events are failing 
the cuts imposed on the reconstructed masses (Eq. (44)). To understand this effect, consider that 
increasing the cm. energy also enlarges the phase space. Consequently, our internal propagators 
are more likely to be on-shell. 

The contribution from the irreducible background for our i^i^jj signal, 

pp -> W^W^W^, pp -> W ± W ± jj, pp -> tt (45) 

events and 

pp jjZZ, pp jjZW, (46) 

wherein leptons from the Z boson escape from a detector, are estimated [24] to be at most a = 
0.08 fb using a comparable list of selection cuts. However, this previous analysis does not impose 
any restriction on the invariant mass of the system as done in Eq. (44), and therefore, realistically, 
the background will be much less than 0.08 fb. In either case, our W' R signal is clearly above 
background. Using a = 0.08 fb as an estimation for our background, we calculate the significance 
and reachability of our W' R signal at the 14 TeV LHC as shown in Fig. 10. With 100 fb^ 1 integrated 
luminosity, a W R signal via the lepton-number violating process can be observed at a 5<r level up to 
a mass of 3 TeV. As evident, the required integrated luminosity for a discovery at the LHC grows 
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Figure 10: Integrated luminosity needed at 14 TeV LHC for (a) achievable statistical significance for 
W' R with Mw = 3 TeV and niM = 500 GeV, and (b) reachable W' R mass at 3cr and 5c sensitivity. 

rapidly with increasing This is expected if we again consider that the W boson becomes 

increasingly boosted as My^i grows. A more boosted W leads to more collimated jets, which have 

R 

more difficulty passing the isolation cuts (Eq. (37)) than their less collimated counterparts. 

V W Chiral Couplings From Angular Correlations at the LHC 

Once a new gauge boson W' is observed at the LHC, it is of fundamental importance to determined 
the nature of its coupling to the SM fermions. Here, we identify various kinematical quantities 
that depend on the chiral couplings of the fermions to a W. Each quantity will have a different 
dependence on the W' chiral couplings and so will provide independent measurements of the chiral 
couplings. 



Figure 11 shows the spin correlations for the process qq' — > W' — > N£ + in the partonic cm. frame for 
both the (a) left-handed and (b) right-handed cases. Double arrowed lines represent spin and single 
arrowed lines momentum. As it is well-known, although the preferred charged lepton momentum 
direction leads to a clear distribution of parity violation, it cannot reveal more detailed nature of 
the chiral coupling. On the other hand, the nature of the W leptonic chiral couplings is encoded 
in polarization of the heavy neutrino, i.e., in the W' R (W' L ) case the heavy neutrino is preferentially 
right-handed (left-handed). Hence, if the polarization of the neutrino can be determined, the left- 
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Figure 11: Spin correlations for qq' — > W' — > Ni + for (a) left-handed and (b) right-handed cou- 
plings. Single arrow lines represent momentum directions and double arrow lines represent spin 
directions. 

handed and right-handed cases can be distinguished. Spin observables such as (sjv ■ «), where sn 
is the spin of the heavy neutrino and a is an arbitrary spin quantization axis, are sensitive to the 
polarization of the heavy neutrino. Defining the angle 6* between the a and the direction of motion 
of the charged lepton originating from the heavy neutrino decay, p£ 2 , the angular distribution of 
the partial width of the neutrino decaying into a charged lepton and two jets is [38] 

:(N -+ e±jj) = \ (l + 2 A l± cos 0*) , (47) 



r dcos 8* 

where A i+ = —A 1 = A due to the CP invariance. The coefficient A is related to (sn • &) and 
is the forward-backward asymmetry of the charged lepton with respect to the direction a. We 
will refer to A as the analyzing power. The angular distribution of either of the two jets from 
the neutrino decay will also have a similar linear form and may be used to perform this analysis, 
although uncertainties in jet measurements may cause more complications. 

A highly boosted neutrino from a heavy W' decay will be produced mostly in a helicity state; 
hence, it is natural to choose a = p^, the direction of motion of the neutrino in the partonic 
cm. frame, and measure p£ 2 in the neutrino rest-frame. At the partonic level, the angular distribu- 
tion of the lepton from neutrino decay in the reconstructable neutrino rest-frame is (See App. B) 



da(ud ->■ ef e t w ~) _ °Tot. 



dcos 9 



/2 



a(W )-a(W T ) \ ( 2-^ N \ / gtf \Y tlN \* - g^V^l 
a(W )+a(W T )J \2 + nl) \ g %\Y ilN \> + &\V hN \ 



Here it(Wo) an d &(Wt) are the partonic level ud — > W' + — > ifi^^x cross sections with N decaying 
into longitudinally (A = 0) and transversely (A = T) polarized W's, respectively. They are given 
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by 



ct(Wo) = a(ud->£fN ->£f{£Wo) (49) 

1 1 g 2 \V^?\VhN? to 2 g 2\ / /2i v |2,^2| F /W\ 
" 9 2*^ (l + ^ 2 ) ^ +5l A^ 1 ^ 1 +^1^1 J^J 



(1 - y^(l - /i^(2 + fijf) -a" ( 5 °) 



[(a - M^,) 2 + (rw/'M^) 2 ] V2y 



ct(Wt) = -> -> ^s^r) ( 51 ) 

= a(W )x2y 2 v . (52) 



where //jv = tun/Vs, yw = Mw/m,N, and &Tot. = (£(Wo) + ct(Wt)) x BR(IF — >• g(f) is the total 
partonic cross section. As W' comes on-shell, /ijy ~~ * X N- I n this reference frame, 6* from Eq. (47) 
satisfies 

cos 6* = cos #£ 2 = j5^ 2 • pn, (53) 

where, again, p£ 2 is measured in the neutrino rest-frame and pn is measured in the partonic cm. 
frame. 

For an on-shell W, the analyzing power at the partonic and hadronic level are the same. In 
such a case, after comparing Eqs. (47) and (48), we find that the analyzing power is 



.4 



1 fa(Wo)-g(W T )\ (2-x 2 N \ f g R 2 \Y hN \ 2 - gj 2 ^-' 2 



iN\ 



2 \a(W ) + a(W T )J \2 + x%J W R 2 \Y eiN \ 2 + 9 e L 2 \V eiN \ 
1 (l-2y 2 w \ (2-x 2 N \ f g R 2 \Y hN \ 2 - gi 2 \V hN \ 2 



(54) 



2 \l + 2y 2 w ) \2 + x 2 N J \g R 2\ YilN \2 + gi 2 \V hN \ 2 , 

The different signs for the analyzing power between the neutrino decays to the two different 
W polarizations and between the W' L R cases can be understood via the spin correlation in Fig. 7. 
For the W' R case, a heavy neutrino decaying to a longitudinal (transverse) W will have the charged 
lepton preferentially moving with (against) p^- For the W' L case the helicity of the neutrino, and 
therefore the direction of the charged lepton, is reversed. Hence the analyzing power is proportional 
to (&(W ) - a{W T )){g R 2 \Y hN \ 2 - <tf\V tlN ?). 

In the analysis of Fig. 7, the left- and right-chiral neutrinos at the W — > N£ + vertex are 
approximated as the left-handed and right-handed helicity states in the partonic cm. frame. As 
the neutrino becomes more massive relative to the W , the approximation of the chiral basis by the 
helicity basis begins to break down, i.e., the left- (right-) helicity state makes a larger contribution 
to the right- (left-) chiral state. In Eq. (48), this is reflected by the cos0£ 2 (cos#£ for simplicity) 
coefficient 

2 ~ X N _ 2M W , - m N 
2 + x 2 N 2M 2 v ,+m 2 N 

As xn increases, the distribution flattens due to the right-handed (left-handed) neutrino helicity 
state, thereby making a larger contribution to the W' L (W' R ) distributions. 
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Figure 12: The angular distribution of the charged lepton originating from neutrino decay in the 
heavy neutrino rest-frame with respect to the neutrino moving direction in the partonic cm. frame 
at the LHC with Mw, set by Eq. (19). Distribution (a) without smearing or cuts, (b) with 
energy smearing and cuts in Eqs. (35), (37), (43), and (44) , and (c) with all cuts applied to (b) 
except the ARjj cuts in Eq. (37). The solid lines are for the Monte Carlo simulation results and 
in (a) and (c) the dashed lines are for the analytical result in Eq. (48). 
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Figure 12 shows the hadronic level angular distribution of the lepton in the neutrino's rest- 
frame for both W' L and W R at the LHC. The case without smearing or cuts is shown in Fig. 12(a), 
and contains both the analytical results (dashed line) and Monte Carlo simulation (solid line) 
histograms. As can be clearly seen, the analytical and numerical results are in good agreement. 
Figure 12(b) shows the leptonic angular distribution after energy smearing and cuts in Eqs. (35), 
(37), (43), and (44). Notice that there is a small depletion of events for cos On ~ 1 and a large 
depletion when cos 0£ < 0. First, when cos On ~ 1 the charged lepton is moving with and the jets 
against the direction of motion of the neutrino in the partonic cm. frame. Hence, with boost back 
to the partonic cm. frame, the jets are softest at this point and the jet pt cuts in Eq. (35) lead 
to a depletion of event in this region. When cos0£ < 0, the lepton is moving against and the SM 
W is moving with the neutrino's direction of motion. Hence, with the boost back to the partonic 
cm. frame, the W is boosted and its decay products highly collimated. Consequently, the ARjj 
cuts in Eq. (37) lead to a large depletion of events. Figure 12(c) shows lepton angular distribution 
with the same cuts as Fig. 12(b) except the ARjj cuts. For comparison, both the Monte Carlo 
simulation with cuts (solid) and analytical results without cuts (dashed) are shown. It is clear that 
the discriminating power of the lepton angular distribution would increases and the Monte Carlo 
distribution approaches the analytical results if the jet isolation cuts are relaxed. 

The analyzing power in Eq. (54) can additionally be related to the forward backward asymmetry 

= a {cos B t > 0)-<j(cos9 e < 0) , . 

a(cas0t > 0) +a(cos0 e < 0) ' [ ' 

Without cuts or smearing, A = A; and for the values of m^r, Myy' stipulated in Eq. (19), 

f + 0.43, W-W a 
[-0.43, W' = W' L 

The simulated values for the forward backward asymmetry with consecutive cuts are shown in 
Table 2. Again, simulations are in good agreement with the theoretical prediction for the forward 
backward asymmetry for no smearing or cuts. As the cuts become more severe, the simulated and 
theoretical values deviate more, however the W' L and W' R cases can still be distinguished clearly. 
Furthermore, as shown in the final row, if the ARjj cuts in Eq. (37) are relaxed, the discriminating 
power of the asymmetry is greatly increased, and the theory and simulation are in much better 
agreement. 

B W' Chiral Couplings to Initial-State Quarks 

Thus far, we have only presented the results to test the chiral coupling of W' to the final state 
leptons. It is equally important to examine its couplings to the initial state quarks. Define an 
azimuthal angle 

^ PN x pg 2 p N x p q 
cos <3> = — — - • — ztt, (58) 

\PNXPe 2 \ \PNXPq\ 
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A 


8 TeV 


14 TeV 


Wr 

L 


it 


Wr 

L 


it 


Reco. without Cuts or Smearing 


-0.42 


0.42 


-0.43 


0.43 


+ Smearing + Fiducial + Kinematics (Eq. (35)) 


-0.46 


0.33 


-0.47 


0.34 


+ Isolation (Eq. (37)) 


-0.11 


0.59 


0.083 


0.72 


+$T + m jj Requirements (Eq. (43)) 


-0.078 


0.62 


0.11 


0.75 


+ Mass Reco. (Eq. (44)) 


0.16 


0.77 


0.18 


0.77 




-0.34 


0.49 


-0.34 


0.49 



Table 2: Forward-backward asymmetry for pp — > W' L R — > fi + ^ + qq' with consecutive cuts at 8 and 
14 TeV LHC. The last row has the same cuts applied as the previous row with the removal of the 
ARjj cuts in Eq. (37). 



as the angle between the qq' — > N£f production plane and N — > W~ decay plane in the neutrino 
rest-frame, where pi 2 is the three momentum of £2, the charged lepton identified as originating from 
the neutrino; pm is the direction of motion of the neutrino in the partonic cm. frame; and p q is the 
initial-state quark momentum. The definition of <3> is invariant under boosts along pjy, hence the 
quark and charged lepton momenta can be evaluated either in the partonic cm. or the neutrino 
rest-frame. The angular distribution between the two planes is thus calculated to be 



da Grot. 



d§ 2vr 



1 + % w ( «W0-*ffi | 

16 2 + u 2 N \a[Wo) + a(W T 




(59) 



The distribution for W' L is 180° out of phase with the W R distribution and the slope only depends on 
the W' chiral coupling to the initial-state quarks. Hence, the phase of this distribution determines 
the chirality of the initial-state quarks couplings to the W' independently of the leptonic chiral 
couplings to the W. 

To understand the distribution in Eq. (59), we consider the spin correlations between the initial 
and final states. As noted previously, the angle $ is invariant under the boosts along p^. So for 
simplicity, we consider the spin correlations in the heavy neutrino rest-frame. Figure 13 shows the 
spin correlations of the neutrino production in the neutrino's rest-frame for both the (a) W' L and 
(b) W' R cases. Like before, single arrowed lines represent momentum directions and double arrowed 
lines spin in the helicity basis. Also, we define the production plane to be oriented in the y — z 
plane such that the y-component of the quark momentum always points along the positive y-axis 
and that z = p^. With this axis convention, $ = —(j)£ 2 , where cj)£ 2 is the azimuthal angle of £2 as 
measured from the positive y— axis. 

Figure 14 shows the spin correlations for the heavy neutrino production and decay with the 
spin quantization axis chosen to be the y direction as defined above. The W' L case is shown in 
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Figure 13: Spin correlations for neutrino production in the neutrino rest-frame. Single arrowed 
lines represent momentum and double arrowed lines represent spin in the helicity basis. The z-axis 
is defined to be the neutrino's direction of motion in the partonic cm. frame and the y-axis is 
defined such that y-component of the initial-state quark momentum is always positive. 

Figs. 14(a,c) and the W' R case in (b,d). The solid dots next to the N and l\ indicate that they have 
no momentum in the y-direction. In the W' R case, the initial-state quark must be right-handed 
and the initial-state antiquark left-handed. Hence, the total spin of the initial-state points in the 
positive y-direction, causing the spin of the neutrino to also point in the positive y— direction. When 
the neutrino decays to a longitudinal or transverse W, the lepton from the neutrino decay has spin 
along or against the y-axis, respectively. For the W' R case, figures 14(b) and (d) show the decay 
into longitudinal and transverse Ws, respectively. Therefore, for the decay into Wq (Wt) case, 
the lepton prefers to move in the same (opposite) direction as the initial-state quark and <I> peaks 
at (±7r). In the W' L case, the direction of motion of £2 relative to the direction of motion of the 
initial-state quark is reversed and the peaks in the $ distribution are shifted by ir. This explains 
the 180° phase difference in the angular distribution, Eq. (59), between the W' L and W' R cases, 
and between the neutrino decay to Wo and Wt- Also, notice that this argument only relies on the 
W' — q — q' coupling and not the W' — N — t chiral couplings. Hence, measuring the distribution 
of the angle between the qq' — > Nl\ production and the N — > l2 + W~ decay planes can determine 
the chiral couplings of a W' to light quarks independently from the chiral couplings of the W' to 
leptons. 

Most of the angular definition and analysis depend on the initial state quark momentum direc- 
tion. Since the LHC is a symmetric pp machine, this is not known a priori. However, at the LHC 
u and d quarks are valence and antiquarks are sea. Hence, the initial-state quark generally has 
a larger momentum fraction than the initial-state antiquark; and the initial-state quark direction 
can be identified as the direction of motion of the fully reconstructed partonic cm. frame. Similar 
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Figure 14: Spin correlations in the neutrino rest-frame as described in Fig. 13. Double arrowed 
lines represent spin with y being the quantization axis and single arrowed lines are the y component 
of the particles. 



techniques have been used for studying forward-backward asymmetries associated with new heavy 
gauge bosons [14,39]. 

Figure 15 shows the $ distributions at the 14 TeV LHC with M' w = 3 TeV for both W' L and W' R . 
From Eq. (59), the amplitude of the $ distribution depends on the ratio tun /Myy, and therefore 
increase m^r to 1.5 TeV. The solid line is the $ distribution with the initial state quark moving 
direction identified as the partonic cm. frame boost direction; the dashed lines is the theoretical 
distribution given in Eq. (59); and in (a) the dash-dot lines are the Monte Carlo truth, i.e. using 
the known direction of the initial-state quark. 

Figure 15(a) does not include cuts or smearing; as can be seen, the Monte Carlo truth and 
theoretical calculation agree very well. The reconstructed distribution has a smaller amplitude 
than the theoretical distribution due to the direction of the initial-state quark being misidentified. 
Figure 15(b) shows the theoretical prediction and reconstructed distribution with smearing and 
the cuts in Eqs. (35,37,43,44) applied. For $ = 0, the SM W is maximally boosted and its decay 
products are maximally collimated. Consequently, the ARjj cut in Eq. (37) causes a large depletion 
of events in the central region. Figure 15(c) shows the reconstructed distribution with the same 
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Figure 15: $ distributions at the 14 TeV LHC with M' w = 3 TeV and m N = 1.5 TeV for fully 
reconstructed events (solid), the analytical result in Eq. (59) (dashed), and Monte Carlo truth 
(dash-dot). Figure (a) is without energy smearing or cuts, (b) with energy smearing and cuts in 
Eqs. (35), (37), (43), and (44), and (c) with the same cuts as (b) without the ARjj cut in Eq. (37). 



cuts as (b) minus the ARjj cut. With the relaxation of this cut, the W' L and W' R cases become 
reasonably discernible with the W' L distribution nearly the same as the theoretical prediction. The 
continued depletion of events at <I> = and $ = ±7r are due to the rapidity cuts on leptons and 
jets, respectively. 



VI Unlike-Sign Dilepton Angular Distributions 

Intrinsically, Majorana neutrinos can decay to positively or negatively charged leptons, and there- 
fore also contribute to the L-conserving process 

PP -> w -> efe^jj. (60) 

These events can be reconstructed similarly to the method described in Sec. IV. However, the SM 
backgrounds for this process, particularly pp — > Zjj, will be larger. Our purpose here is not to 
do a full signal versus backgrounds study, but to comment on the differences between the like-sign 
and unlike-sign lepton cases. Again, ud has a larger parton luminosity than du, so we focus only 
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on W ,+ production: 

pp -»• W' + -»• AT^+ -»• jj (61) 

A W' Chiral Coupling from Angular Distributions 

For the unlike-sign case, we mimic our entire like-sign analysis and reconstruct the polar angular 
distribution of the lepton originating from neutrino decay in the heavy neutrino rest-frame (App. B). 
Respectively, the polar and azimuthal distributions are similar to those in Eqs. (48) and (59) up to 
a opposite sign in front of the angular dependence. 



da a Tot. 



d cos 0£ 2 2 

da &Tot 



a(W ) + a(W T )J \2 + ^ N J W R 2 \Y tlN \* + <f L 2 \V tl 



iN 



d<& 2vr 



_ 3tt 2 m f a(W ) - a{W T ) . , a 
16 2 + »% \a{W ) + a{W T ) 1 




(62) 



(63) 



To understand why the sign of the slope for the L-conserving distributions differ from the 
L- violating distributions, we turn to spin correlations. For W' + , the spin correlations for ud — > 
W' + —> N£ + are shown in Fig. 11 without yet specifying iV's decay. However, we only need to 
analyze the angular correlation in the neutrino decay. The spin correlations are simply obtained by 
replacing the right-handed antilepton in Fig. 7 with a left-handed lepton. Since the direction of the 
spin of the lepton is completely determined by the neutrino spin, which is unchanged between the 
two cases, the effect of the helicity flip is to reverse the direction of the final state lepton momentum 
relative to the z direction. Therefore, the slopes of the lepton angular distribution are opposite for 
the like-sign and unlike-sign lepton cases. These same arguments can be made to show that the 
phases of the $ distribution in Eqs. (59) and (63) differ by 180°. 

The analysis of the two cases also reveals that, unlike the angular distributions, the total 
cross section is independent of having like-sign or unlike-sign leptons in the final state. This may 
be understood by recognizing that the difference between the two final states is tantamount to 
a charge conjugation. Having integrated out the angular dependence, the total cross section is 
invariant under parity inversion. Consequently, by CP-invariance, the total rate is invariant under 
charge conjugation. This behavior is evident in Eq. (29) and Fig. 2, which show that N decays to 
i + W~ and i~W + equally. 

VII Summary 

The nature of the neutrino mass remains one of most profound puzzles in particle physics. The 
possibility of its being Majorana-like is an extremely interesting aspect since it may have far- 
reaching consequences in particle physics, nuclear physics and cosmology. 
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Figure 16: For the opposite sign lepton case, the angular distribution of the charged lepton origi- 
nating from neutrino decay in the heavy neutrino rest-frame with respect to the neutrino moving 
direction in the partonic cm. frame at the LHC with Mw, set by Eq. (19). Distribution (a) 
without smearing or cuts, (b) with energy smearing and cuts in Eqs. (35), (37), (43), and (44) , 
and (c) with all cuts applied to (b) except the ARjj cuts in Eq. (37). The solid lines are for the 
Monte Carlo simulation results and in (a) and (c) the dashed lines are for the analytical result in 
Eq. (48). 
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Given the outstanding performance of the LHC, we are motivated to study the observability 
for a heavy Majorana neutrino TV" along with a new charged gauge boson W' at the LHC. We first 
parameterized their couplings in a model-independent approach in Sec. II and presented the current 
constraints on the mass and coupling parameters. 

We studied the production and decay of W and iV at the LHC, and optimized the observability 
of the like-sign dilepton signal over the SM backgrounds. We emphasized the complementarity of 
these two particles by exploiting the characteristic kinematical distributions resulting from spin- 
correlations to unambiguously determine their properties. Our phenomenological results can be 
summarized as follows. 

1. The heavy neutrino is likely to have a large R.H. component and thus the W' R would likely 
yield a larger signal rate than that for W' L , governed by the mixing parameters as discussed 
in Sec. II. Under these assumptions, we found that at the 14 TeV LHC a 5a signal, via the 
clean channels l^i jj, may be reached for M w > = 3 TeV (4 TeV) with 90 fb _1 (1 ab _1 ) 
integrated luminosity, as seen in Fig. 10. 

2. The chiral coupling of W' to the leptons can be inferred by the polar angle distribution of the 
leptons in the reconstructed neutrino frame, as seen in Fig. 12, owing to the spin correlation 
from the intermediate state N. 

3. The chiral coupling of W to the initial state quarks can be inferred by the azimuthal angular 
distribution of the neutrino production and decay planes, as seen in Fig. 15. 

4. The kinematical distributions for the like-sign and unlike-sign cases have been found to be 
quite sensitive to spin correlations and are complimentary. In particular, the angular distri- 
butions differ by a minus sign and provide qualitative differences for a Majorana and a Dirac 
N . Thus in addition to observing final states that violate lepton-number, comparison of the 
two scenarios provides a means to differentiate the Majorana nature of N. 

Overall, if the LHC serves as a discovery machine for a new gauge boson W', then its properties 
and much rich physics will await to be explored. Perhaps a Majorana nature of a heavy neutrino 
may be first established associated with W' physics. 
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A Neutrino Mixing Formalism and W Couplings 
Al. Neutrino flavor mixing 

Without specifying a mass generation mechanism, we assume that there are three left-handed (L.H.) 
neutrinos (denoted by v a L, a = 1)2,3) with Dirac masses (m), and n right-handed (R.H.) neutrinos 
(denoted by N a iR, a' = 1, . . . , n) with Majorana masses {m'). The mixing between chiral states 
and mass eigenstates may then be parameterized [24] by 



V L \ ( ^3x3 Vs 

N C L ) ~ { X nx3 Y nxn J I N c m , 



(64) 



— T 

where tjj c = Cip denotes the charge conjugate of the spinor field ifj, with C labeling the charge 
conjugation operator, and the chiral states satisfy ipl = (i^ c )l = (^r) - Expanding the L.H. and 
R.H. chiral states, we obtain: 

3 n+3 3 n+3 

m=l m'=4 m=l m'=4 

3 n+3 3 n+3 

V aR = E ^ U ma + E N m' V m'a, N a'R = E ~^ X ™' + E N ™' Y m'a>- (66) 

m=l m'=4 m=l m'=4 

Under this formalism, one expects diagonal mixing of order 1, 

UrfandYY* ~ 0(1); (67) 

and suppressed off-diagonal mixing, 

VV ] and XX ] ~ 0{m/rri). (68) 

A2. Model-Independent W' Charged Current Couplings 

The goal of this paper is to explore the feasibility of quantifying the properties of a new charged 
gauge boson, W', at the LHC. For this purpose, we relax the W interactions to include both 
left-handed and right-handed leptons, 

*■-(?).• (69) 

with a,b = 1, 2, 3. The L.H. neutrinos and charged leptons that are members of SU(2)i doublets 
in the Standard Model (SM) are denoted by v a i and l a . The R.H. neutrinos, which are SM 
singlets, and R.H. charged leptons are denoted by N^r and If,. To grasp the feature of Left-Right 
symmetric models for a W', we pair N^r and lb into the an SU(2)/j doublet. Though there may 
be more "sterile" neutrinos, i.e., b > 3, we consider only 6 = 3 and one new mass eigenstate in our 
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phenomenological presentation. The mass mixing matrix in Eq. (64), in the present case, becomes 
a 6 x 6 matrix with several repeating entires. 

With this assignment, the resulting charged current interactions are 



V2 



(70) 



a=l v 6=1 / 

We have explicitly included the couplings of left- and right-charged currents with new gauge inter- 
actions via W' L R . 

The gauge state leptons, l a and lb, may be rotated into the mass eigenstates, which are defined 
to be the flavors eigenstates I = e,fj,, r. This amounts to the rotation 



(71) 



l=e 



With the SM-like simplest Higgs mechanism, this transformation is trivial and we will make it 
implicit without loss of generality. By simultaneously expanding into the neutrinos' mass basis and 
into the charged leptons' flavor basis, we obtain 

3 n+3 

91 



i=e 



Y.^ U *mi+ E NC m> V r 
m=l m'=4 
3 n+3 

K^mi + E 



m'i 



,m=l 



m'=4 



yw + h.c. 
-fp R r + h.c, 



(72) 



where 



U m£ = E U maO a i, Vrf t = ^ V^, a O ae , X ml = ^ X* lb O b l, Y* t , t = ^ Y m 'bObl- (73) 
o=l a=l 6=1 6=1 

These are the general couplings for the W' charged currents that we follow in this study. 
Leptonic couplings to the SM W boson can be recovered from Eq. (72) by identifying W >± ->■ W ± 
and by setting 

g{ = g and g e R = 0, (74) 

where g is the SU(2)^ coupling constant in the SM. Similarly, we arrive at the SU(2)# charged 
current coupling by identifying W' — > W R and by setting 

g{ = and g R ^ 0. (75) 

In the quark sector, we do not plan to go through a fully- fledged construction for the charged 
current couplings. Instead, we take the simplest approach and just parameterize the model- 
independent W' Lagrangian by 



C = "7= E K + ^ V S KM '^ [9l P L + 9rPr\ di + h.c. 
V m=i 



where y CKM ' is an unknown flavor mixing matrix. 
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Figure 17: The partonic- level process for heavy W + production and decay into like-sign leptons 
and quarks in hadronic collisions. The longer, black arrow not touching the Feynman diagram 
denotes fermion flow (FF). 

B Derivation of Partonic Level Angular Distributions 

We strive clarify a few subtleties that arise when calculating observables involving Major ana 
fermions. To do so, we present a detailed derivation of the matrix element for the lepton-number 
(L) violating process: 

Ui{p A ) + Mpb) -> W'+ -> £f( Pl ) + £+(p 2 ) + q m (p 3 ) + q n (p 4 ), (77) 

with an intermediate Majorana neutrino of mass m^, and governed by the Lagrangian given in 
Section II. As discussed in Section IV, and shown in Fig. 3, there are two interfering Feynman 
diagrams associated with our 2£ + 2j final state. The interference term may be neglect safely when 
calculating the amplitude squared, |A^| 2 , since the heavy neutrino's width is very narrow and 
thus the interference is expected to be small. When constructing and evaluating |.M| 2 , we focus 
on only a single diagram (Fig. 17) but stress that the two diagrams can be treated identically. 
Additionally, the narrowness of the SM W boson's width allows us to further apply the Narrow 
Width Approximation (NWA). The NWA stipulates that, due to its small width compared to its 
mass, the W boson will dominantly be produced on-shell, and further implies 

afadj -> tt l 2<ri) ~ a i u idj -> £i4 W ~) x BR(W -> qq'), (78) 

where BR(X — > Y) is the branching fraction of X going into Y . Since BR{W — > qq') is well-known, 
our work is reduced to determining the analytical expression for 

a(uidj ->£fe£w-). (79) 
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Bl. Determination of the Spin- Summed, Polarization-Dependent, 

Squared Matrix Element 

The usefulness of Feynman rules stems from the ability to assign specific multiplicative factors to 
each component of a Feynman diagram. However, Dirac field Feynman rules are dependent on 
Wick's Theorem, which is a statement on field contractions. For Dirac fields, only combinations of 
the form ipijj can contract, where as for Majorana fields, ipip and tptfj are allowed to contract. In 
short, Feynman rules for Dirac fermions do not account for all possible Majorana interactions. 

We therefore adopt the Feynman rules developed in Ref. [36] for a two-fold reason. The first is 
that the rules for diagram segments not involving Majorana fermions do not change. The second 
is that for parts that do involve Majorana fermions, the new Feynman rules reduce to (a) treating 
the Majorana fermion like a Dirac fermion and modifying the vertex factor for an ordinary Dirac 
fermion with an appropriately placed factor of —1, and/or (b) making a single u -H- v spinor 
substitution. The placement of the additional minus sign and possible spinor substitution is based 
on the direction of fermion flow (FF) relative to the traditionally chosen fermion number flow 
(FNF). When the fermion flow and fermion number flow are equal, the newer rules simplify to the 
usual rules. Computationally, these rules provide a desirable technique that can be automated in 
a straight forward manner. 

In the present case, we identify the relevant FF as being identical to the lepton number-changing 
current. The FF current starts at £±, the charged lepton produced in the W boson decay, and points 
anti-parallel to ^i's momentum; the current then continues parallel to the Majorana neutrino's 
momentum; and finally terminates at £2, the charged lepton produced in the iV decay, and points 
parallel to momentum. See the curved black arrow in Fig. 17. With this orientation, the FF is 
parallel to the FNF at the W'(.\N vertex, and anti-parallel to it at the Ni^W vertex. This change 
in relative current orientation causes two modifications, the first of which is to the spinor of the 
outgoing lepton originating from the Nl^W vertex: 



and accounts explicitly for the change in lepton number. The second modification is to the Nl^W 
vertex itself and occurs in the following way: 



where g is the SM SU(2)l coupling constant, Pr } l = |(1 ± 75), and, as defined in Ref. [36], the 
primed-vertex convention indicates 



(80) 




(81) 



r' = cr T c~ 



1 



(82) 
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where C is the charge conjugation operator and for which 



n 



-l, re{7V'" / } 



(83) 



As a result, we find that the matrix element describing the Uidj — > if^^x scattering process, 
for an outgoing SM W~ boson with polarization A, and in the Feynman Gauge, is 



iM\ = e* Xp {p w ) 



VBjA^UAi ■ u 2 C'fi e2 (fl N + m N )B fJ ,e 1 



(s - M$y, + iTw'Mw') (p% -m 2 N + iT N m N ) 



where the vertex terms are given by 



A% = ^ KM 'r [9 q R ?R + 9lPL] 

B v llN = ±=-f [g e R P R Y £lN + glP L v; iN 



(84) 

(85) 
(86) 



To be explicit: £*\ p {pw) denotes the outgoing polarization vector of the on-shell W boson with 
momentum p\y, mass Myy, and polarization A; v R j represents the the spinor v of an initial-state 
antiquark of flavor j and momentum p R ] similarly, UAi represents the spinor u of an initial-state 
quark of flavor i and momentum pa', u 2 denotes the spinor of our final-state antilepton with flavor 
£2 and momentum p 2 ; and likewise, v\ denotes the spinor of our final-state antilepton with flavor £± 
and momentum p\. The W' mass, width, and momentum-squared are respectively given by Mw, 
T\y/, and the Mandelstam variable 



(87) 



S = (PA+Pb) 2 = (Pi +P2+Pwf- 

The heavy neutrino's mass, width, and momentum are similarly given by uin, Tjy, and 

PN = PA + PB ~ Pi = PW + P2- 

After squaring and summing over external spins, diagrams, and colors (Nc), but not external 
boson polarizations (A), the polarization-dependant squared amplitude is 



4A£ g 2 \vg KM '\ 2 \Vz 2 N? TV 


P'aT 1>Bl^ 


f / R 2 P R + gl 2 P L ) 






2 3 (1 + <W 


{s-M 2 ,,) 2 + {T W ,M W ,) 2 




(p 2 N - m 2 N ) 2 + (r^mAr) 2 



xTr ^ l7(T (p'n + m N ) fa i>2 f x p R (&N + m N ) 7m [g R 2 P R \Y hN \ 2 + g e L 2 P L \V hN \ 



(89) 



9 3/\r2 n 2 \ V CKM'\2 It/ |2 

2 iV c 9 \Vji I \vg 2 N\ 



(1 + 



M^y + {T W ,M W ,) 2 (p 2 N -m 2 N y + {T N m N f 



yhN\ z [g q R g R ) A x + \Y hN \- z [g q L g R ) Bx + \V tlN F[gl<f L ) c x + \v hN { z K<?£ ) ^ 



,1 J 



(90) 
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where 



A x = 2{pa- Pi)(pb • Pn)[{pn -P2) + 2(pat • ex){e\-p2)\ 

- rn 2 N (p A -pi)[(pB ■ P2) + 2(p B ■ £\)(e\- p 2 )], (91) 
B\ = 2(pb -Pi)(pa -Pn) [{pn -P2) + 2(pat • £\){e\ ■ p 2 )\ 

- rn 2 N (p B -pi)[(pA-P2) + 2(pA-£\)(e\-P2)], (92) 
C\ = m 2 N (p A ■ Pi) [(pb -P2) + 2{pb ■ -P2)] , (93) 
Vx = m 2 N (p B -pi)[(pA-P2) + 2(pA-e x )(e x -p2)], (94) 

and £\ is taken to be real. 

The Majorana neutrino's width, Tjy, is expected to be very small. Therefore, to simplify analytic 
integration, we again apply the Narrow Width Approximation such that 

1 vr 



5(p 2 N -m 2 N ). (95) 



(p 2 N - m 2 N ) 2 + {T N m N ) 2 T N m N 

We are motivated to make this additional approximation to highlight and emphasize the analyz- 
ing power of the angular distributions. Our reported numerical results do not reflect this extra 
stipulation (see Eq. (123)). Consequentially, the squared and summed amplitude becomes 

yW M ** N °? \Y$™1 (96) 

) Us-M^,) 2 + {T WI M W ,^ 

\Yi lN \ 2 (&&) 2 Ax + \Y hN \ 2 (g q L g l R ) 2 B x + \v hN \ 2 (glal) 2 C x + \v hN \ 2 (g q R g l L ) 2 v x 
B2. Phase Space Volume Element 

We calculate the partonic-level cross section using the usual formula, 

c 

Here, the factor of 4iVp comes from averaging over initial-state colors and spins. The factor dPS n 
represents the n-body phase space volume element, 

dPS n (P; Pl ...p n ) = f[ {2 ^ k 2Ek (2vr) 4 5 4 (P- Pl p n ) , (98) 

which can be decomposed using the recursion formula 

dPS n (P;pi, . . . ,p n ) = dPS n -i(P]Pi, . . . ,p n -l,n) X dPS 2 (Pn-l.n, Pn-1 , Pn) X '-, (99) 

Z7T 
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where P = Ylm=iPm an d Pi,j — Pi +Pj- I n the present case, dPSs is expressible as 

d tP 

dPS 3 (pA + PB',Pi,P2,Pw) = dPS 2 (pA+ Pb;pi,Pn) x dPS 2 (pN;P2,Pw) x rT~~~ • (100) 

Since each dPSk is individually Lorentz invariant, the two phase space elements in Eq. (100) 
can be evaluated in different reference frames. When dPS2{pi,PN) is evaluated in the partonic 
cm. frame and dPS2(p2,Pw) i n the neutrino rest-frame, the full volume element is found to be 

dPS 3 ( PA + P b;pi,p2,pw) = W 1 "^ x dn j 1 '^ 1 x (101) 

with 

,2 m 2 



„2 _ m N ~2 _PN 2 _ M W , , 



S 



Pn 



and, in the on-shell limit, 



A^Wj hn x N = — — , p w ~^yw = • (103) 

The solid angle element d£l]\[ is defined as the angle made by N with respect to the direction of 
propagation of the initial-state quark in the cm. frame; dQ,£ 2 is defined as the angle made by l\ 
with respect to the heavy neutrino spin axis in the neutrino's rest-frame. 

B3. Partonic-Level Angular Distributions 

The angular distribution of the charged lepton from the neutrino decay is most efficiently determined 
by evaluating ^ |7W| 2 in the neutrino rest-frame. Like individual dPSk volume elements, \Ai\ 2 is 
separately Lorentz invariant and thus can be evaluated in its own reference frame. 

In order to evaluate Eq. (96) in the neutrino rest-frame, we must first rotate and boost the 
four-momenta of the initial-state quarks from the cm. frame. Without the loss of generality, we 
assume that the initial-state (anti)quark is originally traveling in the positive (negative) z— axis 
and that the £fN pair propagate in y — z plane. This allows us to rotate the entire 2 — > 2 system 
such that the neutrino's momentum is aligned with the i— axis, and then boost into the neutrino 
rest-frame. Since we are applying the NWA and immediately integrating over dp N , we will take iV 
to be on-shell. After boosting, our four-momenta are: 

PA = ~r~ — ((1 - cos#tv) + fi 2 N (l + cos6n), 0, -2/iTv sin6>Ar, p? N {\ + cos On) - (1 - cos6>at)) , 

Pb = -. ((1 + cos6>tv) + fi%(l - cos On), 0, 2/ijvsin^jv, A*at(1 — cosOn) — (1 + cos0/v)) , 

4mjv 

p N = (m N , 0,0,0), and Pl = - fi 2 N ){l, 0, 0, -1), (104) 

2m n 
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where On represents the polar angle between and pa in the cm. frame. In the neutrino rest- 
frame, the N — > ltW~ decay products' momenta are 



Pw 



\p2 1 ( 1 , sin 0£ 2 cos (f)£ 2 , sin 0£ 2 sin <f>£ 2 , cos 0£ 2 ) 
E w 



\P2\ 



\P 



m N 



\P2\ 



sm t>£ 2 cos q>£ 2 , — sin U^ 2 sm < 



cos^ 2 I , E\y 



m N 



(l + y 2 w ), (105) 



\P2\ V 2 

where 0£ 2 and <f>£ 2 are defined with respect to the neutrino spin axis in the cm. frame. Explicitly, 
z = pat, where pm = Pn/\pn\ is measured in the cm. frame, and (f>£ 2 w.r.t. to the +y axis. This is 
consistent with Eq. (101). The polarization vectors for the SM W boson are subsequently: 

Ew f\P2\ 



e^ipw) 



■ sm t)£ 2 cos q>£ 2 , — sm t>£ 2 sm <p£ 2 , — cos I 
mw \E\y 



(0, — cos 0£ 2 cos ( 



cos 0£ 2 sm (p£ 2 , sm I 



(0, sin <j) i2 , - cos 4>£ 2 , 0) . 



(106) 



Here the labels 0, Tl, and T2 denote the longitudinal and transverse polarizations of the outgoing 
vector boson. After combining Eqs. (96), (97), (101), and integrating over dp 2 N , as well as dil^, 
for the L-violating process Uidj — > £fN — > t[l\W~ with a longitudinally polarized W~ boson the 
angular distribution is 

a(W ) 



dap 
d£l£ 2 



2 4 



7T 



x {4 



1 + 



o / I 2 q2 

37T fi N 9r~ 9l 



2-/4A ( 9 e R 2 \Y ilN \ 2 - gj 2 \V£ lN \ 2 
2 + ^7 W R 2 \Y£ lN \*+gi 2 \V eiN \ 2 

sm0£ 2 cos <^ 2 }. 



COS 0£ 



Accordingly, for transversely polarized W bosons the angular distributions are 



(107) 



d&Tl d&T2 

dn 



a(W 7 



C-2 



dftf 



2 5 



7T 



x {4 



,,2 \ /„£2|v 1 2 ^2 IT/ 1 2 

/% W #R l^iivl - 5l I^iJvI 



2 + M^; V^l^ivP+^I^ATl 



COS I 



3?r pat / gj; - 



sin 0£ 2 cos <^ 2 }. 



(2 + ^)U 2 + ^f, 
In the preceding lines, we have used the following quantities 



108) 



„2 U/CKM'\2 It/ 1 2 

9 \Vji I \Ve 2 N\ 
3N C 21° ^2 (l + ^ l£2 ) 

771 JV S 



^ 2 + <?rj ^i^r+^iw 

(l-/i Ar ) 2 (l-y^) 2 (2 + ^) 



<l 2 



^ [(I - m 2 ,,) 2 + ovm^) 2 

<t(Wt) = a{ud^£fN^e+l+Wr) 
= a{W )x2y 2 v 



2 V W 



(109) 
(110) 
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Integrating over the azimuthal angle, the polar distributions are calculated to be 



dan 



a(W ) 



d cos ( 



1 + 



N 



2 + » 2 N J \g R 2 \Y £lN \ 2 + g[ 2 \V h 



COS I 



and 



d&T _ d{&Ti + <5ti) <t(WV) 



dcos l 



dcos i 



2-^ \ f <#iW-i£ 2 iv^l 2 

2 + ^7 W R 2 \Y ilN \ 2 + 9 £ L 2 \V ilN \ 2 



cos % 



(111) 



;ii2) 



After combining the two, we find that the polarization-summed polar distribution for the full 



Ujd 



ifi^ll' process is 



da T , 



at. 



dcos 9 



(2 



d(a + a T ) 
d cos 6f n 



<?Tot. 



1 + 



a(W )-&(W T ) f2-fi%\ { g R 2 \Y ilN \ 2 - g e L 2 \V ei 



a(W ) + a(W T ) \2 + ^ N J \g R 2 \Y hN \> + g[ 2 \V ilN 



iN 



where 



a(W ) - a(W T ) _ a{W ) - 2y 2 w a(W ) _ 1 - 2y 2 w 
a(W ) + a(W T ) &(W ) + 2y* r &(W ) l + 2y^ : 
and the total partonic-level cross section is 



cos Of 



(113) 



(114) 



OTot. 



a(uidj — > l\l\qq) 

(ct(Wo) + a(W T )) x BR{W qq) 

a{W Q )(l + 2y 2 w ) xBR(W^qq') 

2 UrCKM'\2 It/ |2 

9 \Vji I \Ve 2 N\ 



3N C 2 1 " ^ (1 + S hi2 ) 



q 2 

9r 



9 q L 2 ){g^\Y, lN \ z + 9V\v h 



N 



m N 



r 



N 



g (1-^(1-^(2 + ^) ^1 + 2^ 



[(« - M^,) 2 + (r^M^) 2 ] V 2^ 



Ml" 



x BR(W qq) 



(115) 
(116) 
(117) 
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Having instead chosen to integrate first over the polar angle before the azimuthal angle, the 
polarization-dependent azimuthal distributions for the Uidj — > £^N — > t^l^W - process are 



dan 



2tt 



1 



3^ 



' <j2 o2' 

9r -9\ 



16 (2 + /4) \g q R 2 + g q L 2 , 



COS ( 



and 



dax _ (daTi +<5"T 2 ) (t(Wt) 



2tt 



1 + 



3tt 2 



' q 2 <7 2 ' 



COS ( 



(119) 



(120) 



16 (2 + ^) \g q R 2 + g q L 2 , 

Similarly, after combining the azimuthal distributions, the total polarization-summed azimuthal 
distribution for the full Uidj — > l^l^qq? process is 



da 



Tot. 



9f, 



2vr 



3vr 2 



a(W ) - a{W T ) 



16 {2 + tf N ) \a(W ) + a(W T )J \g q R 2 + g q L 



' q 2 o2' 

9r -91 



COS i 



(121) 
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Under the definition of the azimuthal angle, <E>, in Eq. (58), we have $ 
recover Eq. (59): 



and consequentially 



dcfTot. _ &Tot. 
d<Z> 2vr 



1 



3vr 2 



fJ-N 



q2 q2^ 



a(W )-a(W T )\ gjC-A 



16 (2 + ^) \a(W ) + a(W T )J \g q R 2 + gf 



cos $ 



(122) 



Lastly, were the NWA never applied to N, the differential cross section for the Uidj 
process is 

da 



g 2 \ V CKM> ? |^|2 



dp 2 



N 



q 2 , q 2\ 2 £2\\r |2, 2 it? \2 I 2 

9r +9l I [Pn9r \ y £iN\ +m N \Ve lN \ g L 



3N C 2 11 vr3 M 2 , (l + 8 hh ) 

s (1 - ^) 2 (2 + A|) (1 ~ Pw) 2 {l + 2^) 

[(J - m^,) 2 + (rvi\%,) 2 ] [Cp& - O 2 + OVm^ 



(123) 



where p? N = p 2 N /s and p\ = M^/p 



B4. Partonic-Level Angular Distributions: L-Conserving Case 

For comparison, we consider the case where the heavy neutrino decays through the following L- 
conserving process: 

-+ W £fN l+qW + . (124) 

Following the identical arguments specified in the preceding appendix, the subsequent polarization- 
dependent angular distributions are 
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(125) 



<if^ 2 dft 
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The polarization-summed distributions for the polar and azimuthal cases are therefore 



(126) 
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respectively, where o~Tot. is still given by Eq.(118). Comparison to Eqs. (113) and (122) demonstrates 
that the slopes of the angular distributions differ in sign for the L- violating and L-conserving cases. 
Consequentially, adding the L— conserving and L— violating distributions together results in the 
quantitative feature 



~ da Tot. da Tot. 

OTot. = -, 7. h 



d cos 9f n d cos 1 



d °Tot. + d °Tot. 



d$ d$ 



(129) 



where L denotes the lepton number-conserving (violating) angular distributions. 
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